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ABSTRACT 

Context. The cosmic infrared background (CIRB) consists mainly of the integrated light of distant galaxies. In the far-infrared the 
current estimates of its surface brightness are based on the measurements of the COBE satellite. Independent confirmation of these 
results is still needed from other instruments. 

Aims. In this paper we derive estimates of the far-infrared CIRB using measurements made with the ISOPHOT instrument aboard the 
ISO satellite. The results are used to seek further confirmation of the CIRB levels that have been derived by various groups using the 
COBE data. 

Methods. We study three regions of very low cirrus emission. The surface brightness observed with the ISOPHOT instrument at 90, 
150, and ISOyt/m is correlated with hydrogen 21 cm line data from the Eff'elsberg radio telescope. Extrapolation to zero hydrogen col- 
umn density gives an estimate for the sum of extragalactic signal plus zodiacal light. The zodiacal light is subtracted using ISOPHOT 
data at shorter wavelengths. Thus, the resulting estimate of the far-infrared CIRB is based on ISO measurements alone. 
Results. In the range 150 to ISOyum, we obtain a CIRB value of 1.08+0.32+0.30 MJy sr~^ quoting statistical and systematic errors 
separately. In the 90fim band, we obtain a 2-cr upper limit of 2.3 MJy sr~^ . 

Conclusions. The estimates derived from ISOPHOT far-infrared maps are consistent with the earlier COBE results. 
Key words, galaxies: evolution - cosmology: observations - infrared: galaxies 



1. Introduction 

The extragalactic background light (EEL) consists of the inte- 
grated light of all galaxies along the line of sight with possi- 
ble additional contributions from intergalactic gas and dust and 
hypothetical decaying relic particles. It plays an important role 
in cosmological studies because most of the gravitational and 
fusion energy released in the universe since the recombination 
epoch is expected to reside in the EEL. Measurements of the 
cosmic infrared background, CIRE, help to address some cen- 
tral, but still largely open astrophysical problems, including the 
early evolution of galaxies, and the entire star formation his- 
tory of the universe. An important issue is the balance between 
the UV-optical-NIR and the far-infrared backgrounds; the frac- 
tion of optical radiation lost by dust obscuration re-appears as 
dust emission at longer wavelengths. The absolute level of the 
CIRE, the fluctuations in the CIRE surface brightness, and the 
resolved bright end of the distribution of galaxies contributing to 
the CIRE all provide strong constraints on the models of galaxy 
evolution through difl'erent epochs. For reviews, see Hauser & 
Dwek (12001 ) and Lagache, Puget, & Dole (l2005l). Th e full anal- 
ysis of the data from the DIREE (Hauser et Si\.l9W, Schlegel et 
al. 1998 ) and FIRAS (Fixen et al. 1997) experiments indicated 



a CIRE at a surprisingly high level of ~1 MJy sr"^ between 140 
and 240 yum. Preliminary results had been obtained by Puget et 
al. ( 1996 ). Lagache et al. ( 1999 ) claimed the detection of a com- 
ponent of Galactic dust emission associated with warm ionised 
medium. The removal of this component led to a CIRE level of 
0.7 MJysr"^ at 140 yum. 
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* Based on observations with the Infrared Space Observatory ISO. 
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UK) and with participation of IS AS and NASA. 



Because the FIR CIRE is important for cosmology these re- 
sults need to be confirmed by independent measurements. The 
ISOPHOT instrument (Lemke at el. I1996K flown on the cryo- 
genic, actively cooled ISO satellite, provided the capabilities for 
this. The ISOPHOT observation technique was difl'erent from 
COEE: (1) with its relatively small f.o.v. ISOPHOT was capa- 
ble of looking into the darkest spots between the cirrus clouds; 
(2) ISOPHOT had high sensitivity in the important FIR window 
at 120-200 yum; (3) with its good spatial and multi- wavelength 
FIR spectral sampling ISOPHOT gave an improved possibility 
of separating and eliminating the emission of Galactic cirrus. 
The primary goal of the ISOPHOT EEL project is the determi- 
nation of the absolute level of the FIR CIRE. The other goals are 
the measurement of the spatial CIRE fluctuations and the detec- 
tion of the bright end of the FIR point source distribution. The 
bright end of the galaxy population contributing to the FIR CIRE 
signal was analysed by Juvela et al. (120000 . 
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2. The method 

We examine three regions of low cirrus emission that were 
mapped with the ISOPHOT at 90, 150, and ISOyum. Because of 
the high sensitivity of the ISOPHOT FIR detectors, we can di- 
rectly correlate HI with ISOPHOT measurements for each FIR 
band separately. In the case of DIRBE, the original analysis per- 
formed by the DIRBE team used lOOyum as an ISM template 
and, therefore, the accuracy of the CIRB detections at 140yum 
and 240 yum also depended on the systematic uncertainties of the 
100 yum data (Hauser et al.f T998l Arendt et al. fT998b . 

The HI lines are optically thin and their intensity traces the 
amount of neutral hydrogen along the line-of- sight. The level of 
FIR emission associated with the ionised medium is still uncer- 
tain and we will consider the possible effects later in the analysis. 

As a first step, a relation between the HI line area and the 
FIR surface brightness is obtained. The relation depends on the 
gas-to-dust ratio, grain properties, and the radiation field illumi- 
nating the interstellar medium (ISM) along the line-of- sight. No 
significant variations have been observed in the gas-to-dust ratio 
apart from those associated with large scale metallicity varia- 
tions. Similarly, because of the diff'use nature of the HI clouds, 
no small scale changes in the intrinsic dust properties or dust 
temperature are expected. Under these conditions the FIR sig- 
nal should have a linear dependence on the HI column density. 
Because each field is considered individually, possible diff'er- 
ences in the HI-FIR relation towards diff'erent regions can and 
will be taken into account. For each field, an extrapolation to 
zero HI intensity eliminates emission associated with the neutral 
ISM (for details, see Sect. 14.11) . The remaining signal is equal 
to the sum of the zodiacal light (ZL) and the CIRB. These com- 
ponents are not removed because they are uncorrelated with the 
HI emission. Furthermore, the ZL has a smooth distribution and 
remains practically constant within each of the areas covered by 
individual ISOPHOT maps (see Abraham et al. .l997) . If the ZL 
level is known, the absolute value of the CIRB can be obtained. 
The ZL estimation is described in detail in Sect. 14. 2[ 



3. Observations 

We study three low surface brightness fields that are labelled 
NGP, EBL22, and EBL26. The field NGP is located at the North 
Galactic Pole, the field EBL22 is similarly at a high ecliptic lati- 
tude, while the third one, EBL26, lies close to the ecliptic plane 
(see Table IB. 3b . EBL26 was selected as a field with high ZL 
level with the purpose of estimating the ZL contribution at the 
diff'erent wavelengths observed in this project. 

The observations of the hydrogen 21 cm line were made 
with the Eff'elsberg radio telescope in May 2002. The telescope 
beam has a FWHM of 9 arcminutes. The areas mapped with the 
ISOPHOT instrument were covered with pointings at steps of 
FWHM/2. The stray radiation was removed with a program de- 
veloped by R Kalberla (see Kalberla [T982l Hartmann et al. T9961 
Kalberla etal. 12005]). 

For details of the observations of the EBL fields and the 
associated data reduction, see Appendix |B1 The principles of 
ISOPHOT data reduction and calibration of surface brightness 
measurements are explained in Appendix lAl 



Table 1. Parameters of linear fits of FIR surface brightness ver- 
sus the HI line area. The l-cr error estimates determined with 
the bootstrap method are given in parentheses. For NGP, the re- 
sults correspond to a fit to the combined data of the northern and 
southern sub-fields (see Appendix IB. 21) . 



Field 


A 


Offset 


Slope 




(jim) 


(MJy sr-^ 


(10-3 MJy sr-i K'^ km'^ s) 


EBL22 


90 


5.53 (0.30) 


35.15 (4.26) 


EBL22 


150 


3.56 (0.34) 


38.91 (4.38) 


EBL22 


180 


3.10(0.38) 


27.95 (5.42) 


EBL261 


90 


18.66 (1.22) 


17.28 (6.57) 


EBL261 


150 


6.38 (1.33) 


27.42 (7.36) 


EBL261 


180 


6.00(1.30) 


25.76 (7.41) 


NGP 


90 


6.31 (0.17) 


24.60 (3.98) 


NGP 


150 


3.53 (0.23) 


28.61 (5.58) 


NGP 


180 


3.20 (0.19) 


31.91 (4.22) 



ipit to data with W(HI) <200 Kkm s"^ only. 



4. Analysis and results 

4.1. Subtraction of Galactic cirrus emission using HI data 

The FIR surface brightness was correlated at each observed 
wavelength with the integrated line area of the HI spectra. At 
each observed HI position the average FIR signal was calculated 
using spatial weighting with a gaussian with FWHM equal to 9\ 
Only those pointings are used where the centre of the Eff'elsberg 
beam falls inside the FIR map. In addition to the observational 
uncertainties, each data point was weighted in direct proportion 
to the fraction of the HI FWHM beam that was covered by FIR 
observations. Therefore, the data close to FIR map boundaries 
get lower weight in the following analysis. 

The obtained correlations are shown in Fig. \T\ For FIR ob- 
servations the plotted error bars are based on the statistical un- 
certainties reported by the PIA. 

The figures include linear fits that take into account the esti- 
mated uncertainties in both FIR and HI data. The slopes and zero 
points of the fit are given in Table [T] In field EBL26 there is a 
clear break in the relation above W(HI)=200Kkms"^ that may 
indicate the presence of molecular gas. There is also one fairly 
bright galaxy that is located in the region of higher cirrus emis- 
sion and may have aff'ected the correlation. Therefore, in the field 
EBL26 the linear fitting was carried out using only data below 
W(HI) = 200 K km s"^ . In the other fields the hydrogen column 
densities are in general smaller, W(HI) ^ 100Kkms"^ so that 
the fraction of molecular gas can be expected to be insignificant. 

The off'sets thus obtained correspond to an extrapolation to 
zero HI column density. To the extent to which the remaining 
contributions of ionised and molecular gas can be ignored (see 
below), the values correspond to the sum of CIRB and the zodi- 
acal light. 

4.2. Subtraction of the zodiacal light 

The zodiacal light (ZL) emission is assumed to have a pure black 
body spectrum. The colour temperature of the spectrum depends 
on the ecliptic coordinates of the source and the solar elongation 
at the time of the observations. Leinert et al. (2002} have studied 
the variations of mid-infrared ZL spectra over the sky using a set 
of observations made with the ISOPHOT spectrometer. We use 
their results to fix the colour temperature of the ZL spectra. 

The absolute intensity of the ZL emission in the FIR is es- 
timated with the help of shorter wavelength ISOPHOT observa- 
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Fig. 1. FIR surface brightness as a function of HI line area W(HI) in the three EBL fields, EBL22 (left), EBL26 (middle), and NGP 
(right). Each point corresponds to one pointing of the HI observations. The uncertainties in the HI line area are estimated based on 
the noise in velocity channels outside detected HI emission. For each HI spectrum the corresponding average FIR signal has been 
calculated using for weighting a gaussian with FWHM=9'. The corresponding error bars are based on error estimates reported by 
PIA from which the formal uncertainties of the weighted mean are calculated. The long dashed line shows the result of a linear fit 
that takes into account the uncertainties in both variables. The dotted lines indicate 67% confidence intervals that are obtained with 
the bootstrap method. 



tions made using the ISOPHOT P detector in the absolute pho- 
tometry observing mode PHT-05 (Laureijs et al. 2003 ). Because 
the observations were made in regions of low cirrus emission, 
the mid-infrared signal is completely dominated by the ZL. The 
measurements were carried out close to the larger raster maps, in 
terms of both time and position. Therefore, they give a good es- 
timate for the zodiacal light emission present in the raster maps. 
FIR absolute photometry measurements were made at the same 
time and at the same positions. These are used to make a correc- 
tion for the contribution that the interstellar dust has, conversely, 
on the measured mid-infrared values. The complete list of obser- 
vations is given in Table IB. 2l 

The derived ZL values obtained from the fits (ZL+cirrus) 
are listed in Table O The values are given at the nominal wave- 
lengths assuming a spectrum vly =constant. The uncertainties 
were estimated based on the quality of the fits (see AppendixiDl). 
In fields EBL26 and NGP, because error estimate of each of the 
two measurements is itself uncertain, we conservatively take the 
average of the two error estimates as the uncertainty of the mean. 



Table 2. The estimated zodiacal light emission. The columns 
are: (1) name of the EBL field (see Appendix, Table IB. 21) , (2) 
temperature of the zodiacal light spectrum (Leinert et al. 2002|), 
(3)-(6) estimated intensity of the zodiacal light at 25 yum, 90yum, 
150 yum, and 1 80 yum, and (7) relative uncertainty of the ZL value 
calculated on the basis of the diff'erence of the fitted ZL model 
and the observations in the range 7.3-90yum. The zodiacal light 
estimates are given at the nominal wavelengths assuming a spec- 
trum y/y ^constant. For EBL26 and NGP, two separate positions 
were observed (see Fig. IB. H and Table IB. 2l) . 



Field 


TzL 


7(25) 


7(90) 


7(150) 


7(180) 


rms 




(K) 




(MJy/sr) 




error 


EBL22 


280 


40.51 


7.69 


3.03 


2.19 


13% 


EBL26_ZL1 


270 


106.05 


21.03 


8.34 


6.02 


21% 


EBL26_ZL2 


270 


100.50 


19.84 


7.87 


5.68 


7% 


EBL26 (aver.) 


270 


100.98 


19.94 


7.91 


5.71 


14% 


NGP_ZL1 


260 


28.65 


5.95 


2.38 


1.72 


5% 


NGP_ZL2 


260 


29.15 


6.06 


2.42 


1.75 


18% 


NGP (aver.) 


260 


28.69 


5.96 


2.38 


1.72 


12% 



4.3. Estimated CIRB levels and their uncertainties 



Table [3] lists the CIRB levels that are estimated based on the 
linear fits between FIR and HI data (Table [B and the zodiacal 
light values of Tabled The uncertainties are obtained by adding 
in quadrature the estimated errors of the off'sets from Table [T] 
the errors of the zodiacal light values from Table [31 and the error 
resulting from the dark current subtraction (see Appendix O, 



^ ^offset + ^L + C^DC- 



(1) 



The uncertainty due to the dark current is estimated to be 
ctdc =0.25-0.30 MJy sr"^ and it is likely to be the main factor 



aff'ecting the uncertainty of the zero point of the FIR observa- 
tions (see Appendix [C])- 

The results obtained for the three individual fields can be 
combined, deriving our final estimates for the CIRB and its un- 
certainty. In the case of the field EBL26 the values are rela- 
tively unprecise because of the high ZL level. This uncertainty 
is reflected in the error estimates. Combining the results we 
get average values -0.54+0.65 MJy sr"\ 0.83+0.41 MJy sr"\ 
1.26+0.37 MJy sr'^ at 90yum, 150yum, and ISOyum, respectively, 
as given in the last line of Tabled 
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Table 3. Estimated level of the CIRB for the individual fields. 
The values correspond to the diff'erence between the off'sets listed 
in Table [T] and the ZL estimates of Table [2l No colour correc- 
tion was applied. The l-cr estimated statistical errors are given 
in parentheses (see text). 



Field 


I(90jim) 


/(150//m) 


/(180yum) 




(MJy sr-^ 


(MJy sr-^ 


(MJy sr-^ 


EBL22 


-2.16(1.04) 


0.53 (0.52) 


0.91 (0.47) 


EBL261 


-1.28 (3.05) 


-1.53 (1.73) 


0.29 (1.53) 


NGP 


0.35 (0.74) 


1.15 (0.37) 


1.48 (0.28) 


average 


-0.54 (0.65) 


0.83 (0.41) 


1.26 (0.37) 



^Fit to data with W(HI) <200 Kkm s"^ only. 



The 90 yum values are very low, because in both the EBL22 
and EBL26 fields negative values are obtained. In the case of 
EBL26 the negative value is not surprising, because the expected 
CIRB level is only a small fraction of the zodiacal light which 
itself has a considerable statistical uncertainty. Therefore, the re- 
sult is sensitive also to any systematic errors of the ZL estimates. 
Apart from the results at 90 yum, the variation between fields is 
only slightly larger than expected on the basis of the quoted er- 
ror estimates. At 150yum a negative value is obtained for EBL26 
which, nevertheless, is less than 2-cr below the highest values. 

At 90yum we can derive for EBL only an upper limit. The 
ISOyum and ISOyum bands are close to each other and the CIRB 
values should be very similar. Therefore, based on the three 
fields and the two frequency bands, we can calculate, as a 
weighted average, an estimate for the CIRB in the range 150- 
ISOyum. The result is 1.08 ±0.32 MJy sr'^ The result would not 
change significantly (less than l-cr) even if either EBL22 or 
EBL26 were omitted from the analysis. 

4.4. The reliability of the CIRB values 

In addition to the statistical uncertainties, the results are aff'ected 
by systematic errors. The CIRB estimates are not aff'ected by 
the HI antenna temperature scale. However, the presence of un- 
subtracted stray radiation could aff'ect the HI zero point of the 
HI data and, thus, lower the CIRB values. We cannot directly 
estimate the presence of residual stray radiation in the HI data. 
However, in Sect. IB. 3 1 we compare some of our HI spectra with 
data from the Leiden/Dwingeloo survey (Hartmann & Burton 
[T9971 Kalberla et al. 120051) and we find that the residual stray 
radiation is likely to be less than 4Kkms"^ which, assuming 
a slope of 29x10"^ MJy sr"^ (Kkms"^)"^ (see Table[T]), corre- 
sponds to 0.12MJysr"^ In this case HI stray radiation would 
not be a major source of error. 

The relative calibration of the ISOPHOT FIR cameras and 
the P-detectors directly aff'ects the estimated FIR ZL levels and 
is probably the most important source of systematic errors. 
According to the ISOPHOT Handbook (Laureijs et al. I2QQ3]) 
the absolute accuracy of the CI 00 and C200 cameras and the 
P-detectors is typically of the order of 20%. If there were a dif- 
ference of 10% in the relative calibration of the MIR and FIR 
bands, this would cause a similar percentual error to the ZL esti- 
mates. The fact that we obtained negative CIRB values at 90 yum, 
especially when the absolute level of the ZL is high, suggests that 
the FIR ZL levels may have been overestimated. The eff'ect of 
an error of 10% would range from ~2MJy sr"^ at 90yum in the 
field EBL26 to ~0.2MJysr-i at 150-180yum in the field NGP 
Taking into account the relative weighting of the three fields, a 
10% error in ZL corresponds to an uncertainty of 0.3 MJy sr"^ 



in the 150-180yum CIRB estimate. Assuming a systematic un- 
certainty of this magnitude, the CIRB estimate can be written as 
1.08+0. 32+0. 30MJy sr"^ where the first error estimate refers to 
statistical and the second to systematic uncertainties. 

At 90 yum the negative value obtained for EBL26 carries very 
little weight and an additional 10% systematic uncertainty in the 
ZL would correspond to an additional uncertainty of 2 MJy sr"^ . 
In EBL22 the CIRB values was -2.16+1.04 and a 10% system- 
atic error in the ZL values would correspond to 0.77 MJy sr"^ 
The CIRB value is ~2-(T below zero and suggests that the ZL val- 
ues may contain a systematic error of -10-20% that has reduced 
the obtained CIRB values. In the field NGP the 90 yum CIRB es- 
timate was 0.35+0.74 MJy sr"^ Assuming a 10% systematic un- 
certainty in the ZL and adding the error estimates in quadrature, 
the CIRB estimate becomes 0.35 +0.95 MJy sr~^ and we obtain 
a 2-(T upper limit of 2.3 MJy sr"^ . 

5. Discussion 

5.1. Dust emission associated with ionised gas 

Our analysis was based on the correlation of HI emission and 
FIR intensity. So far we have omitted the possible eff'ect that 
dust mixed with ionised gas might have. The ionised component 
can aff'ect the results only as far as it is uncorrected with the HI 
emission. 

Lagache et al. (l2QQQb decomposed the DIRBE FIR inten- 
sity into components correlated with the neutral and the ionised 
medium. The column density of ionised hydrogen, N(H'^), was 
estimated based on the Ha line. They found that the infrared 
emissivity of dust associated with the ionised medium would be 
very similar to the emissivity of dust within the neutral medium. 
However, Odegard et al. (2007 ) recently re-examined this issue, 
obtaining significantly lower emissivity values for the ionised 
medium. The derived 2-cr upper limit for the lOOyum emissivity 
per hydrogen ion was typically only 40% of the emissivity in the 
neutral atomic medium. 

We use the all- sky Ha map produced by Finkbeiner (2 0031) 
to examine the possible contribution of the ionised medium to 
the FIR emission. The resolution of the Ha data is 6' for fields 
EBL22 and EBL26 and one degree at the location of the field 
NGP. The average Ha emission in the EBL22, EBL26, and NGP 
fields is -0.7 R, -0.5 R, and -0.6 R, respectively. The Ha back- 
ground contains small scale structure that may be caused by faint 
point sources, mainly stars. Therefore, the quoted Ha levels are 
not caused by the diff'use ISM only. For example, in NGP the 
Ha image is dominated by an unresolved (- one degree) emis- 
sion peak at the centre of the field, the nature of which remains 
unknown. Apart from this, the Ha background does not show 
any significant gradients or correlation with the FIR emission. 
Therefore, we consider only the eff'ect on the average FIR sig- 
nal. Using the Lagache et al. (2000) conversion factors an Ha 
signal of 0.6 R would correspond to -0.5 MJy sr"^ in the FIR. 
Therefore, the CIRB values could be overestimated by a similar 
amount. However, adopting the Odegard et al. (120071) l-cr upper 
limits, the contribution from the ionised medium should remain 
below O.lMJy/sr. Furthermore, in our analysis we have corre- 
lated the FIR emission only with HI while in the quoted studies 
the FIR signal was correlated simultaneously with both HI and 
H^. Therefore, since HI and H^ are themselves correlated, the 
correction to be applied to our results should be correspondingly 
smaller. Therefore, we believe that the possible eff'ects due to 
the presence of an ionised medium are small compared with the 
other uncertainties given above. 
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5.2. Dust emission associated with molecular gas 

If molecular gas is present, the HI lines will underestimate the 
total column density of gas and, because the fraction of molecu- 
lar gas increases with column density, the relation between FIR 
emission and the HI intensity becomes steeper. Our fields have 
low column density and, therefore, the fraction of molecular gas 
should be low. Hydrogen molecules cannot survive in clouds 
with visual extinction below Ay ~ 0.1"^ and, consequently, no 
molecular gas should exist in clouds with column densities be- 
low N(H) ~ 2 X lO^^cm-^. Arendt et al. ([T9981) detected a 
steepening in the FIR vs. HI relation which, however, in diff'er- 
ent regions took place at diff'erent column densities. The eff'ect 
could start already around N(H) ~ 2 x 10^^ cm"^ which corre- 
sponds to an HI line area of W(HI) ~ lOOKkms"^ Kiss et al. 
(12003b observed a change in the spatial power spectra of FIR 
surface brightness around N(H) ~ 10^^ cm"^. This was simi- 
larly interpreted as a sign of the transition between atomic and 
molecular phases. In the EBL fields, molecular emission could 
be significant only in the EBL26 region, where the slope between 
HI and FIR data appears to change at W(HI) ~ 200Kkms"^ 
(see Fig. [T]). Below this limit there is a good, linear correla- 
tion between the FIR surface brightness and HI line area which 
also shows that toward those positions the fraction of molecu- 
lar gas is low. In the EBL estimation only data below W(HI) = 
200 Kkm s"^ were used. 



5.3. Comparison with earlier results 

The earlier CIRB results in the FIR range are all based on mea- 
surements of the COBE satellite. We have derived our CIRB 
estimates using the ISOPHOT measurements, without relying 
on the COBE data even in the determination of the ZL levels. 
Therefore, our result is the first completely independent CIRB 
estimate after the COBE detections. Table IE. II in Appendix [El 
lists the existing CIRB estimates in the FIR wavelength range. 

In the range 150-1 80yum our value is consistent with the 
DIRBE results at 140 yum. According to Kiss et al. (2006) the 
COBE/DIRBE and ISOPHOT FIR surface brightness values 
agree to within ~15% and, therefore, the diff'erences in the sur- 
face brightness scales are likely to be smaller than our statisti- 
cal uncertainty. In the ZL subtraction, the relative calibration of 
the ISOPHOT-P detector and the FIR cameras could introduce a 
systematic error that has a magnitude comparable to that of the 
statistical uncertainty. The low, even negative CIRB estimates 
obtained at 90 yum suggest that this systematic error causes the 
ZL estimates to be -10% too large. Taking into account our sta- 
tistical and systematic uncertainties at 150-1 80yum, we cannot 
exclude even the highest DIRBE estimates close to 1.5 MJy 

At 90yum our 2-cr upper limit of 2.3 MJy sr"^ is consistent 
with the existing DIRBE results. 

Based on the above values, the galaxies resolved with ISO 
FIR observations account for some 5% of the total CIRB (e.g., 
Juvela et al. 2000; Heraudeau et al. 2004; Lagache & Dole 2001 ; 
Kawara et al. 2004). A stacking analysis of Spitzer measure- 
ments (Dole et al. 2006 ) showed that galaxies detected at 24 jim 
contribute some 0.7 MJy sr"^ to the 160yum sky surface bright- 
ness. Therefore, the results from galaxy counts and measure- 
ments of the absolute level of CIRB are converging, and proba- 
bly more than half of the sources responsible for the CIRB have 
already been identified. 



6. Conclusions 

For the ISOPHOT EBL project far-infrared raster maps were ob- 
tained in selected low-cirrus regions. We have analysed these 
observations and, by correlating the FIR surface brightness with 
HI line areas measured with the Eff'elsberg radio telescope, we 
derive estimates for the cosmic infrared background in the wave- 
length range 90-180yum. We determined the level of ZL emis- 
sion using shorter wavelength ISOPHOT observations, without 
relying on a model of the spatial distribution of the ZL emission 
on the sky. Therefore, our results are independent of the existing 
COBE results. 

Based on this study we conclude the following: 

- At 90yum we derived a 2-cr upper limit of 2.3 MJy sr"^ for 
the CIRB. 

- In the range 150-180yum we obtained a CIRB value of 
1.08+0. 32+0. 30MJy sr"\ where we quote separately the es- 
timated statistical and systematic uncertainties. 

- The accuracy of the results was determined mostly by the 
accuracy of the zodiacal light estimates and the dark signal 
subtraction. 

- Assuming the latest emissivity values of dust associated 
with the ionised medium, the uncertainty related to the pres- 
ence of ionised medium was small compared with the other 
sources of uncertainty. 

Acknowledgements. We thank the anonymous referee for valuable comments. 
This work was supported by the Academy of Finland grants no. 1 15056, 107701, 
124620, and 119641. ISOPHOT and the Data Centre at MPIA, Heidelberg, were 
funded by the DLR and the Max-Planck-Gesellschaft. We thank P. Kalberla for 
his help in the planning of the HI measurements and for performing the stray 
radiation correction of these observations. 



References 

Abraham R, Leinert Ch., Lemke D., 1997, A&A 328, 702 

Acosta-Pulido, J. A., Gabriel, C, and Castaneda, O., 2000 Experimental 

Astronomy 10, 333 
Arendt R.G., Odegard N., Weiland J.L. et al. 1998, ApJ 508, 74 
del Burgo, C., Heraudeau, P., and Abraham, P., 2002 in: Proceedings of the 

Symposium "Exploiting the ISO Data Archive - Infrared Astronomy in the 

Internet Age", ESA SP-5111, Siguenza, p. 339 
DwekE. 1998, ApJ. 501, 643 

Dole H., Lagache G., Puget J.-L., et al. 2006, A&A 451, 417 
Finkbeiner D.R, Davis M., Schlegel D.J. 2000, ApJ 544, 81 
Finkbeiner D.RA. 2003, ApJS 146, 407 
Fixen, D.J. 1997, ApJ, 490, 482 

Hartmann D., Burton W.B. 1997, Atlas of Galactic Neutral Hydrogen, 

Cambridge University Press, 
Hartmann D., Kalberla RM.W., Burton W.B., Mebold, U. 1996, A&AS 1 19, 1 15 
Hauser, M.G., et al. 1998, ApJ 508, 25 
Hauser M.G., Dwek E. 2001, ARA&A 39, 249 
Heraudeau R, Oliver S., del Burgo C, et al. 2004, MNRAS 354, 924 
Juvela M., Mattila L., Lemke D. 2000, A&A 360, 813 
Kalberla RM.W., Mebold U., Reif K. 1982, A&A 106,190 
Kalberla RM.W., Burton W.B., Hartmann D., et al. 2005, A&A, 440, 775 
Kawara K., Matsuhara H., Okuda H., et al. 2004, A&A 413, 843 
Kelsall et al. 1998, ApJ, 508, 44 

Kessler M.F., MuUer T.G., Leech, K. et al., 2003, The ISO Handbook I, Mission 
& Satellite Overview, SAI/2000-035/Dc, Noordwijk, The Netherlands 

Kiss C, Abraham R, Klaas U., Juvela M., Lemke D. 2001, A&A 379, 1161 

Kiss C, Abraham R, Klaas U., et al. 2003, A&A 399, 177 

Kiss Cs., Abraham P., Laureijs R.J., Moor A., Birkmann S.M. 2006, MNRAS 
373, 1213 

Klaas, U., Abraham, P., Acosta-Pulido, J. A., et al., 2001 in: Proceedings of the 
Conference "The Calibration Legacy of the ISO Mission", ISO Data Centre, 
Villafranca del Castillo, Madrid, p. 19 

Lagache G., Abergel A., Boulanger R, Desert EX., Puget, J.-L. 1999, A&A 344, 
322 

Lagache G., Dole H. 2001. A&A 372, 702 



6 



M. Juvela et aL: Determination of the cosmic far-infrared background level with the ISOPHOT instrument (RN) 



Lagache G., Haffner L.M., Reynolds R.J., Tufte S.L. 2000, A&A 354, 247 

Lagache G., Puget J.-L. 2000, A&A 355, 17 

Lagache G., Puget J.-L., Dole H. 2005, ARA&A 43, 727 

Laureijs R.J., Klaas, U., Richards, PL, & Schulz B., Abraham P, 2003, The 

ISO Handbook, Volume IV: PHT - The Imaging Photo-Polarimeter, ESA 

SP-1262 (http://isowww.estec.esa.n1/users/handbook/ i 
Leinert Ch., Abraham P., Acosta-Pulido L, Lemke D., Siebenmorgen R., 2002, 

A&A 393, 1073 
Lemke D., Klaas U., Abolins L, et al., 1996, A&A 315, L64 
Lemke D., Kranz Th., Klaas U., et al., 2001, Proceedings of the Conference "The 

Calibration Legacy of the ISO Mission", ESA SP-481, p. 219 
Li A., Draine B. 2001, ApJ 554, 778 

Odegard N., Arendt R.G., Dwek E., et al. 2007, ApJ 667, 1 10 
Puget J.-L., Abergel A., Bernard J.-P, Boulanger R, Burton W.B., et al. 1996, 
A&A 308, L5 

Schlegel D.J., Finkbeiner D.P, Davis M. 1998, ApJ 500, 525 
Stark A.A., Gammie C.R, Wilson R.W. et al. 1992, ApJS 79, 77 



M. Juvela et al.: Determination of the cosmic far-infrared background level with the ISOPHOT instrument (RN), Online Material p 1 



Online Material 



M. Juvela et al.: Determination of the cosmic far-infrared background level with the ISOPHOT instrument (RN), Online Material p 2 



Appendix A: The principles of surface brightness 
observations with ISOPHOT: data reduction and 
calibration 

The most detailed description of the ISOPHOT instrument, 
its observing modes (so called Astronomical Observation 
Templates, AOTs) and the corresponding data analysis and cal- 
ibration steps is given in the ISOPHOT Handbook (Laureijs et 
al. I2QQ3 ). In the following we describe recent calibration tech- 
niques which are beyond the scope of the Handbook and which 
are essential for the determination of the EBL surface brightness. 

A.1. Absolute surface brightness calibration of ISOPHOT 
observations 

ISOPHOT was absolutely calibrated against a flux grid of celes- 
tial point source standards consisting of stars, asteroids and plan- 
ets, thus covering a fair fraction of the entire dynamic flux range 
from ^100 mJy up to about 1000 Jy. Each detector aperture/pixel 
was individually calibrated against these standards. Therefore, 
the basic ISOPHOT calibration is in Jy pixel"^. 

In order to derive proper surface brightness values in 
MJysr"\ the solid angles of each detector aperture/pixel must 
be accurately known 

_/;r(o,o) 

- — (A.l) 

^ eff 

with B;i being the surface brightness, the efl'ective solid 
angle of the pixel/aperture, the total flux of a celestial stan- 
dard and /p^sf^^(0, 0) the fraction of the Point Spread Function 
contained in the pixel/aperture (i.e. the convolution of the PSF 
with the aperture response) when being centred at position (0,0). 
Hence, • (0, 0) is the flux per pixel. 

ISOPHOT's efl'ective solid angles have been determined by 
2D- scanning of a point source over the pixel/aperture in fine 
steps dx' and dy' and measuring the resulting intensity at each 
measurement point (x|, jj), the footprint, taking into account a 
non-flat aperture/pixel response: 

i j 

If the peak of the point source was located outside the aper- 
ture by ~l/2 of the aperture size, the S/N of the resulting in- 
tensity dropped so much that at this border the summation was 
complemented by a model of the broad band telescope PSF and 
adding up the corresponding PSF fractions f^^f^(x',y{) out to 

psi 1 J 

ilOarcmin assuming a flat response, but taking into account a 
cut by ISO's pyramidal central mirror feeding the 4 instrument 
beams. An example of such a synthetic footprint is shown in 
Fig. (AH 

The values of the solid angles used in PI A VI 1.3 are listed 
in Tables lAJlandlAl2l 

It should be noted that an absolute surface brightness cali- 
bration is more accurate than an absolute calibration of a com- 
pact source of similar brightness, since no background subtrac- 
tion has to be performed, which introduces an additional un- 
certainty. The accuracies quoted in the ISOPHOT Handbook 
(Laureijs et al. 12003b , Table 9.1 for extended emission take 
COBE/DIRBE photometry as the reference. By not referring to 
COBE/DIRBE photometry, the absolute surface brightness cali- 
bration for ISOPHOT's ClOO and C200 array is as good as that 
for bright compact sources, i.e. better than 15%. 



Table A.l. Efl'ective solid angles for the 3 x 3 pixels of 
ISOPHOT's ClOO array for the 6 filters with central wavelengths 
Ac. 



filter A, QCioopixeii i^i...9 

(fim) ' (10-^ sr) 



C_50 


65 


0.3737 


0.4005 


0.3929 






0.4167 


0.3963 


0.4120 






0.4064 


0.3690 


0.3943 


C_60 


60 


0.3570 


0.3925 


0.3816 






0.4087 


0.3897 


0.4044 






0.3980 


0.3608 


0.3827 


C_70 


80 


0.4252 


0.4253 


0.4276 






0.4418 


0.4177 


0.4367 






0.4326 


0.3940 


0.4299 


C_100 


100 


0.4855 


0.4517 


0.4666 






0.4679 


0.4396 


0.4632 






0.4603 


0.4194 


0.4701 


C_105 


105 


0.5006 


0.4588 


0.4764 






0.4751 


0.4454 


0.4698 






0.4678 


0.4260 


0.4804 


C_90 


90 


0.4577 


0.4403 


0.4492 






0.4568 


0.4304 


0.4519 






0.4484 


0.4086 


0.4520 



Table A.2. Efl'ective solid angles for the 2x2 pixels of 
ISOPHOT's C200 array for the 5 filters with central wavelengths 
Ac. 



filter ^ ^c2oopi.eii i^i...4 

(jim) (10-^ sr) 



C_160 


170 


1.782 


1.940 


1.895 


1.781 


C_200 


200 


1.810 


1.996 


1.988 


1.878 


C_180 


180 


1.792 


1.960 


1.927 


1.815 


C_135 


150 


1.759 


1.898 


1.825 


1.711 


C_120 


120 


1.722 


1.843 


1.708 


1.587 



A.2. New calibration products and strategies for PI A V1 1.3 

For the very sensitive analysis needed for the EBL determina- 
tion and, in particular, an absolute surface brightness calibration 
that is as accurate as possible, a number of calibration upgrades 
and new calibration features have been developed and imple- 
mented in PI A VI 1.3. For the ones which are not described in 
the ISOPHOT Handbook (Laureijs et al. 2003), we provide a 
description and examples for the ClOO detector in the following. 
An overview of the individual calibration steps associated with 
diff'erent instrument components is shown in Fig. lA.2l By appli- 
cation of all these steps, instrumental artefacts are minimized, 
the resulting detector signals are homogenized and a high cali- 
bration reproducibility and accuracy is achieved. 

A.2.1 . Detector responsivity calibration 

The absolute photometric calibration of an individual measure- 
ment is performed via a transfer calibration using the internal 
calibration sources. This measures the actual responsivity of the 
detector and provides the absolute signal-to-flux conversion. It is 
a separate measurement of each observation mode by deflecting 
the chopper mirror to the field of view of the internal calibrator 
(Fine Calibration Source, ECS). The illumination level of the in- 
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Fig.A.l. Synthetic (outer part, i.e. green and blue coloured 
areas, modelled) footprints (convolution of the ISO telescope 
PSF with the pixel aperture response) of the 3 x 3 pixels of 
ISOPHOT's ClOO array for the 60/jm broad band filter. The 
solid angles of each pixel are obtained by integration over the 
footprint area. 
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Fig. A.2. Scheme of the ISOPHOT calibration steps associated 
with the diff'erent instrument components. The meaning of the 
abbreviations is the following: BSL = Bypassing Sky Light cor- 
rection, DS = detector Dark Signal, RL = Ramp Linearisation, 
TC = signal Transient Correction, and RIC = Reset Interval 
Correction. 



ternal calibrator was not fixed but adjusted as much as possible 
to the expected brightness level of the sky. This was achieved by 
selecting an appropriate heating power for the internal source. 
There exists a calibration relation between this heating power 
and the optical power received by each detector pixel which is 
established from measurements on celestial standards. 

Therefore, for reliable and accurate transfer calibrations, the 
following requirements are put on the FCS: 

1) High reproducibility. This was better than 1%, since the 
monitoring of the flux of faint standards was reproducible 



within a few percent, and this uncertainty was dominated by 
the signal noise (Klaas et al.'2001). 
2) A very detailed characterization of the illuminated power de- 
pending on the heating power applied to the source. This is 
illustrated in Fig. lA.3l It involves the following steps: 

a) For each ClOO and C200 array filter all measurements of 
celestial standards done in raster map mode were evalu- 
ated such that for each pixel the background signal was 
properly subtracted and the resulting source signal was 
associated with the celestial standard flux. The ratio of 
the source signal and the simultaneously obtained FCS 
signal gave the illumination power by the FCS for the 
selected heating power. The discrete results were fitted 
and the reliable lower and upper heating power limits 
covered by measurements were identified (Fig. |A3] up- 
per left). The heating power ranges were not identical or 
equally large for each filter (Fig. lA.3l upper right). In gen- 
eral they were shifted to smaller heating power values for 
longer wavelengths. 

b) For fine discrete steps in heating power the inband pow- 
ers were read from the relations and were converted 
to monochromatic surface brightnesses by applying the 
bandpass conversions derived from the relative system 
response profiles, see ISOPHOT Handbook (Laureijs et 
al. 2003), section A.2, and the solid angles of Tables I A. II 
and IA.2I These fluxes were fitted with a modified BE 
curve after appropriate colour correction (Fig. |A3] mid- 
dle left). If for a certain filter the selected heating power 
was outside the reliable limits, the value of this filter was 
excluded from the fit. The fit gave the temperature of the 
FCS for the selected heating power. An additional con- 
straint was that the temperature had to be the same for the 
fit curves of all pixels. ClOO and C200 filter values had 
to be fitted independently because of the diff'erent detec- 
tor areas and hence illumination factors, however, the fits 
were checked for consistent temperatures, because the il- 
luminating FCS was the same for both detectors. 

c) This was achieved for the heating power range from 0.07 
up to 6.5 mW adopting an emissivity of the source oc A}-^^ 
yielding the temperature vs. heating power relation as 
shown in Fig. I A. 3 1 middle right. 

d) By applying this FCS temperature model and the estab- 
lished illumination factors for each pixel it was possible 
to establish homogeneous calibration curves of the in- 
ternal reference source, thus polishing out measurement 
outliers aff'ecting the initial empirical curves. The multi- 
filter approach connecting all curves and not treating 
them individually enabled a large extension and a com- 
mon range for all filters: compare Fig. IA.3I lower centre 
with Fig. lA. 3 1 upper right. 



A.2. 2. Bypassing sky light correction of FCS signal 

As a safety design against single point failures, ISOPHOT was 
not equipped with any cold shutter to suppress straylight when 
performing internal calibration measurements. Therefore, when 
deflecting the chopper onto the illuminated internal calibration 
(FCS) sources, some fraction of the power received on the de- 
tector did not come from the FCS but from sky light bypassing 
along non nominal light paths. Since this depends on the sky 
brightness it is subtracted in the transfer calibration measure- 
ments on celestial standards and hence has to be subtracted for 
any FCS measurement in order to get a reproducible zero point. 
This was achieved by performing a number of measurements on 
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Fig. A.3. Steps in the generation of a homogeneous and most complete calibration of ISOPHOT's long wavelength internal calibra- 
tion sources (FCS). This is illustrated for the central pixel (#5) of the ClOO array camera. Upper left: Measured relation between 
optical power received on the detector and the heating power applied to the internal source. Dots indicate the discrete measurements, 
the solid line is a fit. Upper right: Display of the input curves for all ClOO filters within the reliable heating power range. Middle left: 
For a selected heating power (here: 1.0 mW) monochromatic and colour corrected fluxes of all filters are fitted by a modified BB 
curve. Middle right: By repeating the fits with the same modified BB type for the whole heating power range covered the relation 
between heating power and temperature of the internal source is established. Lower centre: By applying the FCS model the relation 
between optical power and heating power is homogenized and extended to the maximum heating power range covered by at least 
one measurement in any of the ClOO or C200 filters. 



the switched-off', i.e. cold FCS, so that only the bypassing sky 
light contribution was measured. 

The result for one ClOO array pixel is shown in Fig. IA.4I 
which demonstrates a linear dependence of the bypassing sky 
light contribution to the FCS signal on the sky background. This 
correction was established for all ClOO and C200 array pixels. 
The bypassing sky light contribution contains the detector dark 
signal contribution, cf . Sect. IA.2.51 

A.2.3. Effective pixel/aperture solid angles 



These are described in the previous section lAJl and their values 
are compiled in Tables lAJl and lA.2l 



A.2.4. Filter profiles 

The bandpass system responses and the conversion factors from 
inband power to a monochromatic flux, as well as colour correc- 
tion factors are described in the ISOPHOT Handbook (Laureijs 
etal.|2003D. 



A.2.5. Detector dark signal 

The detector dark signals were re-analyzed as described in del 
Burgo (2002). In this analysis special care was given to exclude 
those dark measurements sufl'ering from memory efl'ects by pre- 
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Fig. A.4. Bypassing sky light contribution to the FCS signal de- 
pending on the sky background. 



ceding bright illuminations, thus not representing the true dark 
level. An example of the new results is shown in Fig. lA.SI for the 
central pixel 5 of the ClOO array. A slight orbital dependence is 
visible with an increase of the dark signal towards the beginning 
and the end of the observational window. It can also be noticed 
that there is a scatter of the dark signals at the same orbit po- 
sition and there are occasional large outliers. These are not due 
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to signal determination uncertainties, but are real variations due 
to space weather effects on different revolutions over the ISO 
mission. 
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Fig. A.5. Orbit dependent dark signal determination for the cen- 
tral pixel 5 of ISOPHOT's CI 00 array. Dots represent individ- 
ual measurements obtained during the entire ISO mission, filled 
and open signals identify a different reset interval in the inte- 
gration of the dark signal. The solid line is the fit to the mea- 
surements providing the so-called default dark level. The dotted 
line is the default dark level of an older calibration version used 
before 2001. 



A.2.6. Ramp linearisation 

This was performed as described in the ISOPHOT Handbook 
(Laureijs et al. 12003b . For ISOPHOT's far-infrared detectors two 
types of effects cause non-linearities of the integration ramps: 

1) De-biasing effects of the photoconductors operated with low 
bias caused by feed-back from the integration capacitor. 

2) Non-linearities generated in the cold read-out electronics. 

A.2.7. Signal dependence on reset interval correction 

Despite the ramp linearisation step, signals obtained under con- 
stant illumination, but with different reset intervals show a sys- 
tematic difference, see Fig. IA.6I upper panel. In order to have a 
consistent signal handling of measurements with different reset 
interval settings applied - to optimize the dynamic range of the 
signal - all signals were converted as if they were taken with 
a 1/4 s reset interval. The correction relations were established 
from special calibration measurements applying the full suite of 
reset intervals under constant illumination and this for different 
illumination levels. In this way signal corrections were estab- 
lished for all reset intervals in the range 1/32 s to 8 s (Fig. IA.6I 
middle and lower panel). While previously, as still described in 
the ISOPHOT Handbook (Laureijs et al. [20031) a linear correla- 
tion with offset was used, a re-analysis (del Burgo et al. 2002|) 
yielded non-linear relations as shown in Fig. IA.6I This latter 
analysis also found a bi-modal behaviour for CI 00 array pixels, 
such that the pixels on the main diagonal, #1, 5 and 9, behaved 
differently from the rest of the pixels. For the C200 array all pix- 
els behaved in the same way. 



A.2.8. Signal transient correction 

The ISOPHOT detectors were photoconductors operated under 
low background conditions provided by a cryogenically cooled 
spacecraft. Under these conditions they showed the behaviour 
that the output signal was not instantaneously adjusted to a flux 
change but rather, following an initial jump by a certain frac- 
tion of the flux step, the signal adjusted with some time constant 
to the final level, see e.g. Acosta et al. (2000). In particular the 
ISOPHOT CI 00 detector showed significant transient behaviour. 
This time constant depended on the detector material (doping of 
the semi-conductor and its contacts), the flux step, the direction 
of the flux step (dark to bright versus bright to dark) and the 
illumination history. Attempts had been made to model this be- 
haviour (Acosta et al. 2000), but no unique description could 
be found for the FIR detectors. To overcome this effect at least 
partly the method of transient recognition was implemented in 
the ISOPHOT analysis software as described in the ISOPHOT 
Handbook (Laureijs et al. 12003 ) using the most stable part of 
the measurement for signal determination. Finally, another ap- 
proach was to use a data base of long measurements with signals 
stabilising and to determine the deviation from the end level for 
shorter intermediate times (del Burgo et al. 2002), see Fig. IA.7I 
for an illustration. A measurement time of 128 s was used as ref- 
erence, because 

1) Most calibration measurements in staring mode were per- 
formed with this basic measurement time. 

2) In most cases the signals stabilised within this time. 

For the C200 array pixels the signal transient effect is consid- 
erably smaller and faster and therefore it is sufficient to apply the 
transient recognition as described in the ISOPHOT Handbook 
(Laureijs et al. 2003j . 



Appendix B: Observations and data reduction for 
the EBL fields 

at. ISOPHOT observations 

The following tables give details of the ISOPHOT observations 
used in the paper. Table IB. II lists the raster maps and abso- 
lute photometry measurements that were made at 90, 150, and 
ISOyum. Correspondingly, Table IB. 2l lists observations used for 
the determination of the zodiacal light levels. These include both 
mid-infrared measurements carried out with the ISOPHOT-P 
detector and longer wavelength absolute photometry measure- 
ments carried out with the CI 00 and C200 cameras. 

Each field was mapped in the PHT22 staring raster map 
mode (ISOPHOT Handbook, Laureijs et al. 2003) using filters 
C_90, C_135, and C_180. The corresponding reference wave- 
lengths of the filters are 90yum, 150yum, and ISOyum. The 90yum 
observations were made with the CI 00 detector consisting of 
3x3 pixels, with 43.5''x43.5'' each. The longer wavelength ob- 
servations were made with the C200 detector which has 2x2 
detector pixels, with 89.4''x89.4'' each. The same raster maps 
were used in Juvela et al. (120001) . Table IB. 3] lists the coordinates 
and the sizes of the maps. Additionally, we make use of PHT25 
absolute photometry measurements (see ISOPHOT Handbook, 
Laureijs et al. 2003) made at the same three wavelengths. Two 
positions in NGP, two positions in EBL26, and one position in 
EBL22 were observed in this mode. 
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Table B.2. Observations used for the determination of the zodiacal light emission. The columns are: (1) name of the field, (2)- 
(3) position, (4) wavelength, (5) the ISO identifier number (TDT) of the observation, and (6) time diff'erence between the listed 
observation and the observation of the EBL raster maps of Table IbTT] A time diff'erence is quoted only when the observations were 
not performed within the same day. Observations at wavelengths below 60 jim are made with the ISOPHOT-P detector. 



Field 


RA 


DEC 


A 


AOT 


TDT number 


At 




(J2000) 


(J2000) 


(pim) 






(days) 


EBL22 


2 26 34.4 


-25 53 53 


25 


PHT05 


82101132 


+2 








90 


PHT25 


82101126 


+2 








150 


PHT25 


81901927 










180 


PHT25 


81901529 




EBL26_ZLli 


1 20 3.0 


1 32 30 


7.3 


PHT05 


58600239 










25 


PHT05 


58600240 










60 


PHT25 


58600341 










170 


PHT25 


58600444 




EBL26_ZL2 


1 19 30.2 


1 23 59 


3.6 


PHT05 


58600646 










7.3 


PHT05 


58600647 










25 


PHT05 


58600648 










60 


PHT25 


58600349 










90 


PHT25 


58600342 










90 


PHT25 


58600350 










150 


PHT25 


58600443 










150 


PHT25 


58600451 










170 


PHT25 


58600452 










180 


PHT25 


58600545 










180 


PHT25 


58600553 




NGPJZLl 


13 43 16.3 


40 28 47 


7.3 


PHT05 


55800331 


-4 








25 


PHT05 


55800332 


-4 








60 


PHT25 


55800233 


-4 








90 


PHT25 


55800226 


-4 








90 


PHT25 


55800234 


-4 








150 


PHT25 


56200127 










150 


PHT25 


56200135 










170 


PHT25 


56200136 










180 


PHT25 


56300129 


+ 1 








180 


PHT25 


56300137 


+ 1 


NGPJZL22 


13 43 32.4 


40 33 35 


7.3 


PHT05 


55800123 


-4 








25 


PHT05 


55800124 


-4 








170 


PHT25 


56200128 





^The position is outside the raster maps. 

^The position is outside maps other than the 180 yum map consisting of 15x15 rasters (TDT number 21301803). 

Table B.3. The positions and sizes of the observed fields. Columns are: (1) name of the field, (2)-(3) equatorial coordinates of the 
centre of the field, (4)-(5) galactic coordinates, (6)-(7) ecliptic coordinates, (8) number of raster points, (9) area in square degrees, 
and (10) additional remarks. All areas were observed at 90, 150, and ISOyum. In NGP(N) an additional square map was observed at 
ISOyum only. Details of the individual measurements are listed in Appendix, in Table IB. ll 



Field 


Ra (2000) 


Dec (2000) 


/ 


b 


A. 




Rasters 


Area (□") 


Notes 


EBL22 


02 26 34.5 


-25 53 43 


215.0 


-68.7 


23.8 


-37.9 


32x3 


0.19 






02 32 48.0 


-25 54 06 


215.6 


-67.3 


25.5 


-38.5 


31x1 


0.07 


ID scan into a region of higher cirrus 


EBL26 


01 18 14.5 


01 56 40 


136.5 


-60.2 


18.8 


-5.9 


32x4 


0.27 




NGP(N) 


13 43 53.0 


40 11 35 


86.5 


73.0 


184.2 


46.5 


32x4 


0.27 






13 42 32.2 


40 29 12 


87.9 


73.0 


183.6 


46.6 


15x15 


0.56 


180yL/m only 


NGP(S) 


13 49 43.7 


39 07 30 


81.3 


72.9 


186.4 


46.1 


32x4 


0.27 





B.2. Reduction of EBL field observations 

The ISOPHOT data were processed with PIA (PHT Interactive 
Analysis) program version 11.3. For details of the analysis 
steps, see the ISOPHOT Handbook (Laureijs et al. 2003) and 
Appendix, section|Al For CI 00 a method of signal transient cor- 
rection was introduced in PIA 1 1.3. This procedure was used for 
all CI 00 measurements. Nevertheless, some of the internal cal- 
ibrator (FCS) measurements show residual drifts. In those cases 
we applied transient recognition which removes the initial, un- 
stabilised part of the measurements. The flux density calibration 



was made using the internal calibrator measurements (FCSl) 
performed immediately before and after each map for actual de- 
tector response assessment. The calibration was applied using 
the average response of the two FCS measurements. 

The reduced data contained a few artifacts. These include 
short time scale detector drifts at the beginning of some CI 00 
observations, temporary signal variations caused by cosmic ray 
glitches, and occasional drifting of some detector pixels that may 
also be connected with cosmic ray hits. The time ordered data 
were examined by eye. For rasters and detector pixels afl'ected 
by clear anomalies (glitches or drifting) the corresponding PIA 
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Table B.l. List of ISOPHOT observations of EBL fields carried 
out in the PHT-22 and PHT-25 observation modes. 



Field/ 


A 


RA 


DEC 


Size 


TDT 


AOT 


(yL/m) 


(J2000) 


(J2000) 




number 


EBL22 












PHT22 


90 


2 26 34.4 


-25 53 49 


32x3 


82101111 


PHT22 


150 




— " — 


— " — 


81901910 


PHT22 


180 


— — 


— " — 


— " — 


81901509 


PHT22 


90 


2 32 47.9 


-25 54 6.0 


32x1 


82101115 


PHT22 


180 




-"- 


xl32 


81901513 


PHT22 


150 


-"- 


-"- 


xl 32 


81901914 


PHT25 


150 


2 26 34.4 


-25 53 52.8 


1 


81901927 


PHT25 


180 


2 26 34.4 


-25 53 52.8 


1 


81901529 


PHT25 


90 


2 26 34.4 


-25 53 52.9 


1 


82101126 


EBL26 












PHT22 


90 


1 18 14.5 


1 56 39 


32x4 


58600319 


PHT22 


150 




— " — 


" 


58600418 


PHT22 


180 






" 


58600517 


PHT25 


150 


1 17 39.3 


2 3 9.8 


1 


60302067 


PHT25 


180 


-"- 


-"- 


1 


60302275 


PHT25 


90 


1 19 30.2 


1 23 59.2 


1 


58600350 


PHT25 


90 


-"- 


-"- 


1 


58600342 


PHT25 


150 


1 19 30.2 


1 23 59.2 


1 


58600443 


PHT25 


150 


-"- 


-"- 


1 


58600451 


PHT25 


180 


— " — 


— " — 


1 


58600545 


PHT25 


180 


-"- 


-"- 


1 


58600553 


NGP(N) 












PHT22 


90 


13 43 53.8 


40 11 42 


32x4 


55800211 


PHT22 


150 








56200110 


PHT22 


180 








56300109 


PHT22 


180 


13 42 32.2 


40 29 12 


15x15 


21301803 


PHT25 


90 


13 43 16.3 


40 28 46.5 


1 


55800226 


PHT25 


90 


13 43 16.3 


40 28 46.5 




55800234 


PHT25 


150 


13 43 16.3 


40 28 46.6 




56200127 


PHT25 


150 


13 43 16.3 


40 28 46.5 




56200135 


PHT25 


180 


13 43 16.3 


40 28 46.7 




56300129 


PHT25 


180 


13 43 16.3 


40 28 46.5 




56300137 


NGP(S) 












PHT22 


90 


13 49 44.5 


39 07 37 


32x4 


55800215 


PHT22 


150 








56200114 


PHT22 


180 








56300113 



error estimates were scaled upwards, typically by a factor of a 
few. For each detector pixel the signal values were scaled so that 
their average value over a map became equal to the overall aver- 
age over all detector pixels. The scaling takes into account the al- 
ready manually adjusted error estimates. The flat fielding would 
actually not be necessary, because FIR fluxes are compared only 
with observed HI 21cm lines and, therefore, averaged over areas 
that are large compared with the size of the ISOPHOT rasters. 

Long term detector response drifts are not taken out by a 
simple averaging of the FCS measurements, nor is an initial non- 
linear drift corrected for by linear interpolation between the two 
FCS measurements. Both could introduce an artificial gradient 
in the time ordered data and, because of the systematic scan pat- 
tern, also in the maps. The maps were compared with IRAS data 
in order to see if there were any gradients uncorrelated with the 
IRAS lOOyum signal. The only significant diff'erence was found 
in the C200 observations of the southern NGP field. The gradient 
was removed while keeping the average surface brightness un- 
changed. The correction has little eff'ect on the subsequent anal- 
ysis. Apart from the EBL22 field, all maps contain four detec- 
tor scans that run alternatively in opposite directions along the 
longer map dimension. When data are correlated with the lower 



resolution HI observations, the subsequent scan legs tend to can- 
cel out any long term drifts. 

The raster map observations themselves do not contain any 
direct measurement of the dark current. In such cases one usually 
relies on the orbit dependent "default" dark current estimates in- 
cluded in the PIA. However, absolute photometry PHT-25 mea- 
surements were carried out within a couple of hours before or 
after each raster map. The data reduction was carried out also 
using the dark current and cold FCS values obtained from those 
measurements. In the subsequent analysis, we use maps that are 
averages of those obtained using default dark current values and 
those obtained using PHT-25 dark current measurements. 

When absolute photometry points were inside the mapped 
area they were compared with the surface brightness of the raster 
maps. The maps were re-scaled so that the final surface bright- 
ness corresponds to the average of the original FCS calibrated 
maps and the values given by the absolute photometry measure- 
ments. This causes systematic lowering of the surface brightness 
values of the original maps. For EBL26, NGP(N), and NGP(S) 
the change is typically ~4%, for both CI 00 and C200 observa- 
tions. In the case of EBL22 the correction is larger, some 20%, 
for the C200 detector. 

In the region NGP there are separate northern and southern 
fields that overlap by a few arc minutes. The maps, each contain- 
ing 32x4 raster points, were fitted together using the overlapping 
area, where the final map is at a level equal to the average of 
the northern and the southern maps. The resulting change in the 
surface brightness levels of individual maps was ~5% or less. 
In the north there is yet another 15x15 raster map that was ob- 
served only at ISOyum. Because that measurement includes only 
very short FCS measurements, it was scaled to fit the already 
combined long ISOyum map. This required scaling of the surface 
brightness values by a factor of 1.05. 

The main maps of the field EBL22 cover an area of low 
cirrus emission. There are additional one-directional scans that 
extend to a region of higher surface brightness in the west. In 
the absence of scans in the opposite direction, it is not possible 
to directly determine the presence of detector response drifts. 
However, these observations were reduced using the average of 
the responsivities given by the two FCS measurements and the 
error bars reflect also the difl'erence in the responsivity before 
and after the measurement. Using the overlapping area, the 32 x 1 
raster strips were scaled to the same level with the 32 x 3 raster 
maps. The scalings applied were 0.97, 1.02, and 0.84 at 90yum, 
150yum, and ISOyum, respectively. 

The final FIR errorbars show the uncertainty for the weighted 
means over the Eff'elsberg beam. The noise of each HI spectrum 
was estimated separately using the velocity channels outside the 
line. The uncertainty of the line area was calculated assuming the 
same, uncorrelated noise for the integrated velocity interval. This 
might underestimate the total uncertainty, because it ignores the 
uncertainties in the stray radiation subtraction that do not aff'ect 
the signal in the line wings. However, for a small field the stray 
radiation causes a constant systematic error rather than statistical 
uncertainty and does not aff'ect the weighting of the observations 
when the linear fit is made. 

For selected positions there exist mid-infrared observations 
made with the ISOPHOT P-detectors as well as further abso- 
lute photometry measurements with the CI 00 and C200 cameras 
(see Appendix, Table IB.21) . These observations were performed 
for the purpose of estimating the zodiacal light. The data reduc- 
tion of P-detector data is similar to that of the CI 00 and C200 
cameras, except that also signal linearisation is included. 
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Fig.B.l. The ISOPHOT EBL fields. The three frames show the ISOyum and 90yum ISOPHOT maps. The coordinates correspond 
to the ISOyum maps. The 90yum maps cover the same area but, in the figure, the 90yum maps have been plotted south of the ISOyum 
maps. The small yellow circles indicate the positions observed with the ISOPHOT P-detector for the determination of the zodiacal 
light levels. To indicate the locations of the fields with respect to the galactic and ecliptica planes, the positions are shown on an 
all-sky map that is combined from DIRBE observations between the wavelengths of 12/im and 240 //m. 



B.3. HI measurements 

The observations of the hydrogen 21 cm line were made with 
the Eff'elsberg radio telescope in May 2002. The observed posi- 
tions, 580 in number, are indicated in Fig. IB. 21 The integration 
times were 30 s in EBL22 and 62 s in EBL26. In the field NGP 
the observations were done with 62 s integrations except for the 
northern part where the integration time was 94 s. The average 
noise estimated from the velocity intervals outside the HI line 
is 0.15 K per channel of Ikm/s. This corresponds to a typical 
uncertainty of 1.7 Kkms"^ in the integrated line area. 

For calibration purposes and for precise subtraction of the 
stray radiation, regular observations of the standard region S7 
were made. The stray radiation subtraction is crucial because it 
affects the zero point of the estimated HI column densities. The 
observed fields, NGP in particular, have some of the lowest line- 
of- sight column densities over the whole sky. Under these con- 
ditions the stray radiation received by the telescope side lobes 
becomes a significant fraction of the total signal. The stray ra- 
diation was removed with a program developed by P. Kalberla 
(Kalberla et al. 2005; see Sect. [3]). 

In Fig. IB. 31 we compare our data with spectra from the 
Leiden/Dwingeloo survey (Hartmann & Burton 1997; Kalberla 
et al. 120051) . For this comparison, in order to match the resolution 
of the Leiden/Dwingeloo survey, the Eff'elsberg data were con- 
volved with a gaussian with FWHM equal to 36'. The HI profiles 



agree very well. Part of the diff'erences may be caused by the fact 
that our HI maps do not cover the whole area of the 36' beams. 
Nevertheless, the figure shows that the observations and the stray 
radiation subtraction (see Sect. sect:obs) are consistent with the 
Kalberla ( 2005j results. 



Appendix C: Calibration accuracy 

The error estimates listed in Table [T] are based on the statisti- 
cal uncertainties in the fits between FIR and HI data. The scat- 
ter of data points around the fitted lines is usually larger than 
their estimated uncertainty. This could be a sign of underesti- 
mated measurement uncertainties but is more likely caused by 
true scatter in the relation. If the formal uncertainties of the line 
parameters were estimated based on the error estimates of the 
individual points, the uncertainties could be severely underesti- 
mated. Therefore, instead of relying only on the measurement 
uncertainties, the uncertainty of the fit parameters was estimated 
separately with the bootstrap method so that they reflect the true 
scatter of observed points. The error estimates corresponding to 
a 67% confidence interval are given in Table [T] 

These uncertainties do not include estimates for the system- 
atic errors introduced by the independent calibration of each map 
or the absolute accuracy of the overall ISOPHOT calibration. 
There are both multiplicative and additive sources of uncertainty. 
The former include, for example, uncertainties in the internal 
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calibration source (FCS) measurements (e.g., detector drifts) that 
alter the estimated detector response. The uncertainties that af- 
fect the zero point of the intensity scale are more critical, because 
the CIRB is small compared with the observed signal and can be 
recovered only as the residual after the subtraction of the ZL. 

Table IC.l] lists an assessment of uncertainty that, using data 
in Table [H have been converted into uncertainty of the FIR 
flux at zero HI column density. The quoted values are half of 
the difl'erence of two values obtained in two independent ways. 
Thereby the quoted values are also -l-cr estimates for the un- 
certainty of the average of the two values. 

In Table IC.l I column 4 has been obtained by comparing the 
fine calibration source measurements performed before and af- 
ter each map. The numbers indicate the statistical uncertainty of 
the detector response measurements. The FCS measurements are 
generally very consistent, particularly in the case of the C200 de- 
tector. On the other hand, the eff'ect of the drift aff'ecting the first 
FCS measurement of the one-dimensional strip map of EBL22 
is clearly visible at 90/im. 

The dark signal subtraction is the most important correction 
aff'ecting the zero point of the FIR intensity. Close to each of 
the raster map observations, we have one or two absolute pho- 
tometry observations which include dark signal measurements 
of their own. In PI A, the default dark current calibration is based 
on a larger set (-70) of dark current measurements for which the 
orbit trend has been determined. Therefore, the PI A default dark 
current calibration is less aff'ected by the noise of individual mea- 
surements but may not take into account short time scale varia- 
tions in the detector dark current on a specific orbit. The maps 
were reduced using the default dark current values and the actu- 
ally measured dark current values. In Table IC.ll column 5 shows 
the associated uncertainty in the FIR signal at zero HI column 
density. The observed uncertainty in the dark current values is 
comparable with the variation observed in the systematic analy- 
sis of a large sample of ISOPHOT observations (del Burgo et al. 
120021 see also Fig.lA31). 

When absolute photometry measurements existed within 
mapped areas, those were used to re-scale the surface brightness 
values of the maps (see Sect. IB. 21) . The diff'erence in the abso- 
lute photometry and mapping measurements is used to derive 
the values in column 6 of Table IC.ll The final column reflects 
the difl'erence in the surface brightness in areas where two in- 
dependently calibrated maps overlap. The numbers in columns 
6 and 7 include, of course, dark current and FCS uncertainties 
as one of their components. For the CI 00 observations at 90 jim 
the uncertainty is close to 1 MJy sr"\ i.e., comparable with the 
expected EBL signal. On the other hand, for the C200 detector 
the uncertainty of an individual map is ~0.3 MJy sr"^ Most of 
this is caused by the uncertainty in the dark current values. 

C.1. Straylight radiation 

Stray light may be another instrumental artefact afl'ecting the zero 
level of the FIR surface brightness. By design and operation 
ISO's viewing direction stayed by several tens of degrees away 
from the brightest FIR emitters in the sky, the Sun, the Earth 
and the Moon (Kessler et al. 2003 ). A dedicated straylight pro- 
gram was executed verifying by deep "difl'erential" integrations 
that the uniform straylight level due to these sources was below 
ISOPHOT's detection limit, even under the most unfavourable 
pointing conditions close to the visibility constraints (Lemke et 
al. 2001). 

Specular straylight by the second brightest class of objects, 
the giant planets Jupiter and Saturn, was observed when point- 



ing to within 15' to 1° of the planet, expressing itself as finger- 
like stripes or faint ghost rings (Kessler et al. 2003 , Lemke et al. 
2001 ). The NGP and EBL 22 fields are far away from the eclip- 
tic and can thus not suff'er from this type of straylight. For EBL 
26 we checked the positions of the planets Mars, Jupiter, Saturn, 
Uranus and Neptune at the time of the observations, 1997-06- 
26 and 1997-07-1 1, respectively. Mars, Jupiter, Uranus, Neptune 
were all far off'. Saturn was at a distance of 3.25 degrees, which 
is still more than a factor of 3 off' of any known straylight-critical 
distance. 



Appendix D: Determination of the ZL levels 

The ZL level was estimated by fitting ZL and cirrus templates 
to ISOPHOT observations in the wavelength range from 7.3 yum 
to 200 yum (see Sect. 14.21) . Figure ICTI shows the results of these 
fits. In the field EBL22 we had observations of one position and 
in the fields EBL26 and NGP of two positions (see Table lR2l) . 
For the latter two fields, the figures show the fit to data combined 
from the two positions. 

Table IB. 21 lists the time diff'erence between the listed obser- 
vations and the observations of the raster maps. In the case of 
NGP these are relative to the 150yum observations. The 90//m 
maps were observed four days before and the 1 80 yum one day 
after the 150yum maps. According to the Kelsall et al. (1998) 
ZL model the four day diff'erence causes only -1.5% change 
in the expected ZL. The combined NGP map is almost 1.5 de- 
grees long. In the Kelsall model the diff'erence in the centre po- 
sitions of the southern and northern parts corresponds to about 
1% diff'erence in the ZL. Therefore, we use only one zodiacal 
estimate value for both NGP(N) and NGP(S) and for all obser- 
vations made during the five day interval. 

In the fields EBL26 and NGP, MIR observations exist for 
two separate positions (see Fig. IB. II) . In both fields, the mea- 
surements at these two positions are close to each other, both in 
time and position. Therefore, their ZL values should be identi- 
cal and also the cirrus levels should be very similar. Comparison 
of the fits performed using these independent sets of measure- 
ments gives the first indication of the statistical uncertainty of 
the ZL values. In both fields, the ZL values obtained for the two 
positions agree within 10%. 

The observations are fitted as a sum of ZL and cirrus compo- 
nents. The ZL template is a black body curve at the temperature 
obtained from Leinert et al. (120021) . The cirrus template is based 
on the model by Li & Draine flOOT). Using the ISOPHOT fil- 
ter profiles we calculate for both radiation components, ZL and 
cirrus, and for each filter the in-band power values that can be 
directly compared with the observed values. In the fit we have 
only two free parameters, the intensity of the ZL component and 
the intensity of the cirrus component. The ZL estimates should 
be based mainly on data between lOyum and 60yum where the 
ZL is clearly the dominant component. Therefore, in the fit, the 
weight of the data points in this wavelength range is increased 
by a factor of two. 

The level of the cirrus component is determined mostly by 
the longer wavelength data. In reality, the component corre- 
sponds to the sum of the cirrus and CIRB signals. As long as the 
component is small in the MIR, the ZL estimates are almost in- 
dependent of the exact shape of this template. We confirmed this 
by replacing the Li & Draine (2001) cirrus template by a pure 
CIRB template, using the model curve from Dole et al. ( 2006] 
Fig. 13). The resulting change in the ZL estimates was less than 
one per cent. 
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Table C.l. Assessment of the calibration uncertainty for the ISOPHOT maps. Columns are (1) name of the field, (2) wavelength, 
(3) average surface brightness of the map, (4) diff'erence between calibration measurements performed before and after each map, 
(5) diff'erence between actual dark current measurements and default dark current values, (6) diff'erence between the independently 
calibrated absolute photometry measurements and raster maps, and (7) diff'erence between partially overlapping maps. These un- 
certainties have been converted to correspond to the uncertainties at zero hydrogen column density using the fit parameters listed in 
Tableffl 



Field 


A 


<S > 


A(FCS) 


A(DC) 


A(Abs.) 


A(Join) 




(//m) 


MJy sr~^ 


MJy sr~^ 


MJy sr~^ 


MJy sr~^ 


MJy sr~^ 


EBL22 


90 


9.0 


0.12 


0.43 


-1.15 




EBL22 


150 


5.7 


-0.06 


1.44 


-0.64 




EBL22 


180 


4.5 


-0.38 


0.48 


-0.28 




EBL221 


90 


6.5 


-1.09 


0.22 




-0.17 


EBL221 


150 


3.6 


0.02 


0.30 




0.07 


EBL221 


180 


4.5 


0.04 


1.00 




-0.50 


EBL26 


90 


20.6 


-0.09 


-0.80 


-1.66 




EBL26 


150 


4.3 


-0.05 


-0.28 


-0.09 




EBL26 


180 


3.7 


0.06 


-0.20 


0.11 




NGP(N) 


90 


7.8 


0.49 


-0.19 


-0.50 


0.58 


NGP(N) 


150 


5.1 


0.06 


-0.30 


-0.35 


0.19 


NGP(N) 


180 


4.8 


0.22 


-0.22 


-0.26 


0.39 


NGP(S) 


90 


6.9 


-0.32 


0.26 




-0.49 


NGP(S) 


150 


4.7 


0.04 


-0.30 




-0.18 


NGP(S) 


180 


4.7 


0.05 


-0.31 




-0.31 


^The narrow strip map. 



The actual statistical errors of the ZL values are estimated 
using the standard deviation of the relative errors when observa- 
tions are compared with the fitted ZL curve. The last column of 
Table |2] lists the corresponding error of the mean, calculated us- 
ing data points between 7.3 yum and 90yum. In the case of fields 
NOP and EBL26, the error estimates are calculated from the fits 
where we have combined the data from the two measured posi- 
tions within each field. In all three fields, the obtained relative 
uncertainties are ~10%. In the fields EBL26 and NOP the uncer- 
tainties are also consistent with the diff'erence of the ZL values 
obtained for the two individual positions. The ZL fits are shown 
in Fig.lCT] 

In this paper we have used original ISOPHOT observations 
without applying colour corrections. Therefore, in the fitting pro- 
cedure also the ZL and cirrus templates were converted to corre- 
sponding values using the ISOPHOT filter profiles. However, for 
Fig. IC.ll we have performed colour corrections. The templates 
are plotted by connecting the values at the nominal wavelengths 
by a straight line. The template spectra used in the ZL fitting are 
colour corrected using their respective spectral shapes. In the fig- 
ure, the colour correction of the observed surface brightness val- 
ues is done assuming the blackbody ZL spectrum below 90 yum, 
and a modified black body cirrus spectrum, By(T = 18 K)y^, at 
90yum and longer wavelengths. 

The plots include DIRBE values from the DIRBE weekly 
maps. These correspond to the DIRBE pixel closest to the cen- 
tre of the corresponding ISOPHOT map. Linear interpolation 
was performed between the weeks in order to accurately match 
the solar elongation of the ISOPHOT observations. In addition 
to the DIRBE value that corresponds directly to the ISOPHOT 
observations (solid squares) we plot the DIRBE value for the 
same solar aspect angle and opposite sign of the solar elonga- 
tion. Assuming that the zodiacal dust cloud is symmetric along 
the ecliptic, the two values should be identical. The predictions 
of the ZL model of Kelsall et al. (1998 ) are also plotted. The 
DIRBE values are colour corrected. As in the case of ISOPHOT 
data, colour correction of the observations assumes a blackbody 
ZL spectrum at and below 60 yum, and a modified black body cir- 



rus spectrum, By(T = 18K)y^, at the longer wavelengths, 100, 
140, and 240 yum. 

There is a clear diff'erence in the ISOPHOT and DIRBE sur- 
face brightness scales. The DIRBE values are consistently lower 
by some 20-30%, in the MIR range. In the FIR bands the ex- 
tended cirrus structures combined with the much larger pixel size 
and noise in the DIRBE pixels precludes direct comparison. The 
determination of the CIRB values is not directly aff'ected by a 
possible calibration diff'erence between DIRBE and ISOPHOT 
because, in this paper, we use exclusively ISOPHOT measure- 
ments. 

Systematic uncertainties aff'ecting all ISOPHOT bands have 
only little impact on the derived CIRB values. The relative cali- 
bration accuracy between the FIR cameras and the ISOPHOT-P 
photometer is more important, because the zodiacal light esti- 
mates are based on the latter. When the absolute level of the zo- 
diacal light was estimated we calculated the scatter between the 
SED model and the observations at diff'erent wavelengths (see 
TableO. The scatter was typically -10-20%. The importance of 
this error source depends, of course, on the absolute level of the 
ZL emission. The field EBL26 is located near the ecliptic plane 
and at 90yum the observed signal and the ZL are both of the 
order of 20MJysr"^ Therefore, a relative uncertainty of 10% 
would already correspond to about twice the expected level of 
the CIRB. For EBL22 and especially for NGP the zodiacal fight 
level is much lower so that more meaningful limits can be de- 
rived for the CIRB also at 90 yum. 

The quoted ZL error estimates reflect the uncertainty in the 
determined ZL level in the mid-infrared. If there were a sys- 
tematic difl'erence in the calibration of the mid- and FIR-bands, 
the ZL estimates could be wrong by the corresponding amount. 
Generally the relative calibration accuracy is considered to be 
within 15%. This uncertainty would not necessarily be reflected 
in the quality of the ZL spectrum fits, because a systematic cal- 
ibration error could have been partly compensated by a change 
in the intensity of the cirrus component. 

The ZL spectrum was assumed to be a pure black body with 
the temperature given by Leinert et al. (I2QQ2I) . As far as the mid- 
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Table E.l. Comparison of existing CIRB estimates in the FIR 
range. The error estimates quoted by the authors are shown in 



parenthesis. In our case, 


we include only the statistical uncer- 


tainty. 








A 


ly 


Reference 


Instrument 


(jim) 


(MJy sr-^ 






90 


<2.3 


this paper 


ISO/ISOPHOT 


150/180 


1.08 (0.32)2 


this paper 


ISO/ISOPHOT 


100 




Hauser et al. 1998 


COBE/DIRBE 




0.73 (0.20)1 






100 


0.37 (0.10) 


Dweketal. 1998 


COBE/DIRBE 


100 


0.78 (0.20) 


Lagache et al. 2000 


COBE/DIRBE 


100 


0.83 (0.27) 


Finkbeiner et al. 2000 


COBE/DIRBE 


140 


1.49 (0.33) 


Schlegel et al. 1998 


COBE/DIRBE 


140 


1.17 (0.33) 


Hauser et al. 1998 


COBE/DIRBE 


140 


0.70 (0.28) 


Hauser et al. 1998 


COBE/FIRAS 


140 


0.70 (0.28) 


Lagache et al. 1999 


COBE/DIRBE 


140 


1.12(0.56) 


Lagache et al. 2000 


COBE/DIRBE 


140 


1.17 (0.37) 


Odegard et al. 2007 


COBE/DIRBE 


240 


1.36 (0.16) 


Schlegel et al. 1998 


COBE/DIRBE 


240 


1.12(0.24) 


Hauser et al. 1998 


COBE/DIRBE 


240 


1.04 (0.16) 


Hauser et al. 1998 


COBE/FIRAS 


240 


0.88 (0.16) 


Lagache et al. 1999 


COBE/DIRBE 


240 


0.88 (0.56) 


Lagache et al. 2000 


COBE/DIRBE 


240 


1.04 (0.24) 


Odegard et al. 2007 


COBE/DIRBE 



1 Hauser et al. did not claim detection at lOOyt/m, because the CIRB 
signal failed test for isotropy. 

2 Only the statistical error is quoted. 



infrared points are concerned, a wrong temperature would, at 
some level, be reflected also in our error estimate. However, if 
the ZL spectrum deviated from the assumed shape only in the 
FIR this could again be masked by a change in the fitted cirrus 
component without a corresponding increase in the rms value. 
Therefore, we must explicitly assume that the same ZL temper- 
ature is applicable both at mid-infrared and far-infrared wave- 
lengths. However, because a 5 K change in the ZL temperature 
corresponds to only -2% relative change in the ratio of 150yum 
and 25 jim intensities, this source of uncertainty is unimportant 
compared with the uncertainty in the relative calibrations of the 
different detectors. 

Appendix E: Comparison with DIRBE EBL 
estimates 

This present study represents the first determination of the abso- 
lute level of the FIR EBL that is independent of measurements 
of the CORE DIRBE instrument. In Table [Owe fist FIR EBL 
estimates given in seven publications based on the DIRBE mea- 
surements. Included are also our 2-cr upper limit at 90/jm and 
the EBL estimate for the range 150-180yum. 



C100 
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o r ■ ^ ^^^^ 
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C100 




S_8 (V/s) 

Fig.A.6. Correction of the signal dependence on the selected 
reset interval. Upper panel: Demonstration of the eff'ect, show- 
ing the resulting signal versus the selected reset interval over the 
range from 1/32 s up to 8 s (reset intervals were commanded in 
powers of 2) under constant illumination. Middle panel: Solid 
line: Correction relation for a reset interval of 8 s w.r.t. the ref- 
erence reset interval of 1/4 s for all CI 00 array pixels, except 
the ones on the main diagonal. Dotted line: old linear correlation 
used before the re-analysis. Lower panel: Solid line: Correction 
relation for a reset interval of 8 s w.r.t. the reference reset inter- 
val of 1/4 s for all CI 00 array pixels on the main diagonal (pixels 
#1, 5, and 9). Dotted line: old linear correlation used before the 
re-analysis (same as for middle panel). 
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Fig.A.7. Empirical signal transient correction for ISOPHOT's 
CI 00 array. The left column shows the signal loss for integra- 
tion times of 4, 8, 16, 32, and 64 s (commendable integration 
times of ISOPHOT detectors) with regard to the reference time 
of 128 s. The red line is a fit through the measured points over the 
covered signal range and is used as the correction relation. The 
right column shows the residuals after applying this correction. 
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Fig.B.2. The figures shows as black rectangles the areas 
mapped with ISOPHOT (90, 150, and 180yum) and as circles the 
pointings used in the Eflfelsberg HI observations. The diameter 
of the circles, 9', is equal to the FWHM of the Eff'elsberg beam. 
The frames a, b, and c correspond to regions EBL22, EBL26, 
and NGP. In the case of NGP, the dashed red line indicates the 
area that was mapped at 1 80 jim only. 



M. Juvela et aL: Determination of the cosmic far-infrared background level with the ISOPHOT instrument (RN), Online Material p 13 



6 




''^ -60 -40 -20 20 40 

v [km/s] 

Fig. B.3. Comparison of HI spectra from the Leiden/Dwingeloo 
survey (Kalberla et al. 120051 dashed lines) and our Effelsberg 
data convolved with a beam of 36' (solid lines). The spectra 
correspond to positions at the southern and northern end of the 
NGP map (13M2"^0' +40°30'0'' and 13^52"^0' +38°40'0'0 
and one position in the field EBL26 (1^ 17^ 0' +2° 20' 0''). The 
EBL26 spectra have been scaled by a factor 0.5. 
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Fig. C.l. Fits used to estimate the ZL levels in the three fields 
EBL22, EBL26, and NGP (frames a, b, and c, respectively). The 
red circles are ISOPHOT observations. The lower lines are the 
cirrus (blue solid line) and the ZL (red solid line) templates, 
the uppermost solid green line is their sum. The figures also 
show DIRBE values for the closest DIRBE pixel, read from the 
DIRBE weekly maps. The solid squares correspond to observa- 
tions with the same solar elongation as in the case of ISOPHOT 
observations, the open squares to the other measurement with 
identical absolute value of the solar aspect angle but opposite so- 
lar elongation. For clarity, the latter have been shifted slightly in 
wavelength. The dashed line shows the predictions of the Kelsall 
et al. (jTWSl) ZL model. 



